Comment on "Classical Simulations Including Electron Correlations for Sequential 

Double Ionization" 
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In a recent Letter the authors investigate trajecto- 
ries of a two-electron atomic system which models argon 
submitted to a strong, close to circularly polarized laser 
pulse. In the field-free case, the following classical Hamil- 
tonian is considered: 
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where G 1R 3 is the coordinate vector of the i-th electron 
in a frame centered at the nucleus, Pi are the canonically 
conjugate momenta, and = | | and pi — \pi\. Here 
the Heisenberg core potential Vh is given by 
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and is supposed to incorporate the Heisenberg uncer- 
tainty principle in a classical way 0. Using Hamilto- 
nian (UJ , the authors of Ref . [l[ analyze the double ion- 
ization of argon in order to interpret the experimentally 
measured ionization times of Ref. [H, and in particular 
the anomalous ionization time of the second electron. In 
a very striking way, Hamiltonian ([1]) is able to capture 
a shorter than expected ionization time of the second 
electron in excellent agreement with the experiment. Be- 
low, we show that this excellent agreement is only the 
result of tweaking the parameter a in the potential ([2]) , 
and that this result breaks down as a is varied, in con- 
trast with what the authors claim [l[ "the results of our 
calculations do not depend on this parameter." Their 
result is therefore not structurally stable and does not 
hold as a is made arbitrarily large, which is the regime 
where the potential has been designed (for the constraint 
rp > £ to be satisfied) 0. Commonly used values for a 
are 4 or 5, which is a compromise between the constraint 
and the speed of the integration given the stiffness of the 
problem. In Ref. l|, the choice of one of the param- 
eters (£, a) relies on the second ionization potential of 
the atom [E^ — —1-01 for argon). More precisely, the 
minimum of the one electron Hamiltonian (kinetic en- 
ergy plus Coulomb interaction with the nucleus plus the 
Heisenberg core potential) has to be equal to E 2 - As a 
result, £ is chosen equal to 
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The parameter a is chosen arbitrarily, provided that the 
ground state of the one electron Hamiltonian occurs at 
a configuration where rp — £, which requires that a > 
1.256 (the solution of e x — 1 — 2x = 0). This translates 
into a minimum value for £, approximately equal to 1.190. 




FIG. 1. Mean-value of the time delays between the two suc- 
cessive ionizations in Hamiltonian JTJ) for 1 = 4 PW ■ cm" 2 . 
The other parameters are set as the ones in Ref. 0] for a pulse 
duration of 33 fs. The continuous curve is for 3D calculations 
and the dashed curve corresponds to 2D calculations. 



For a — 2, the one electron energy minimum occurs at 
£ ss 1.259 and not at £ = 1.225 as stated in Ref. [lj. 
However choosing a value of £ given by Eq. Q results in 



an empty ground state energy surface at E g 



-1.59 a.u. 



for the two electron system, and thus a different value of 
£ has to be considered. In the following, we consider the 
value chosen in Ref. [lj, i.e., £ = 1.225. 

In Fig. [TJ we represent the time difference between the 
ionization of the two electrons for the parameters consid- 
ered in Ref. Q (laser wavelength 788 nm, ellipticity 0.77, 
pulse duration 33 fs) and for an intensity of 4 PW-cm~ 2 . 
We see that the time difference is in agreement with the 
measured one (on the order of 10 fs) for a — 2, as in 
Ref. [l[. However it increases as a is made larger, up to 
reaching a value comparable to the one obtained using 
soft-Coulomb potentials (on the order of 30 fs), both for 
the 2D and the 3D calculations. Therefore we argue that 
the quantitative agreement on ionization times is only 
coincidental. The relevance of adding the Heisenberg 
core potential has not been demonstrated. Note that 
a quantitative agreement with the experimental results 
has been obtained in Ref. 0] using soft-Coulomb poten- 
tial, i.e., without invoking the Heisenberg core potential. 
We reached similar conclusions by using soft Coulomb 
potentials with varying effective charge for the electron- 
nucleus interaction and by decreasing the influence of the 
electron-electron correlation. 

In order to obtain qualitative agreement with experi- 
mental results (e.g., the knee shape of the double ioniza- 
tion probability versus intensity or the patterns of the ion 
momentum distribution), it has been shown in the past 
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few decades that the specific choice of the potential does 
not matter. This means that the theoretical or numerical 
interpretation of these experiments is structurally stable. 
However using these models to obtain quantitative agree- 
ment is more subtle since it requires a fine tuning of the 
Hamiltonian model of the atom. This problem is central 
whether one works in the quantum or classical frame- 
work. 

C. Chandre 1 , A. Kamor 1 ' 2 , F. Mauger 1 , and T. Uzer 2 
1 Centre de Physique Theorique, CNRS - Aix-Marscillc 
Universite, Campus de Luminy, Case 907, 13288 Mar- 
seille cedex 09, France 



2 School of Physics, Georgia Institute of Technology, At- 
lanta, Georgia 30332-0430, USA 



[1] Y. Zhou, C. Huang, Q. Liao, and P. Lu, Phys. Rev. Lett. 

109, 053004 (2012). 
[2] C.L. Kirschbaum and L. Wilets, Phys. Rev. A 21, 834 

(1980). 

[3] A.N. Pfeiffer, C. Cirelli, M. Smolarski, R. Dorner, and U. 

Keller, Nature Phys. 7, 428 (2011). 
[4] X. Wang, J. Tian, A.N. Pfeiffer, C. Cirelli, U. Keller and 

J.H. Eberly, preprint (2012). 



